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ABSTRACT 

! 

The  snapping  of  a  simply-supported  shallow  sinusoidal  arch 
under  a  sinusoidally  distributed  step  pressure  load  is  considered. 

In  the  presence  of  velocity-dependent  damping  of  any  nonzero  magni-' 
tude,  it  is  shown  that  there  exists  no  difference  between  static  and 
dynamic  snap  loads  for  arch  rises  above  a  certain  magnitude  and  for 
sufficiently  small  external  disturbances.  Below  the  foregoing  value 
of  arch  rise  snapping  is  governed  entirely  by  symmetric  snap-through. 
The  results  obtained  herein,  when  compared  to  an  existing  analysis  on 
the  subject,  imply  that  a  jump  in  the  critical  snap-through  load  occurs 
at  the  boundary  between  the  damped  and  undamped  systems.  Similar 
results  concerning  cylinders  and  rings  are  mentioned. 
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1.  Introduction 


In  several  recent  theoretical  investigations  concerning  the 
dynamic  stability  of  elastic  shell-type  structures  subjected  to  trans¬ 
ient  loads,  it  has  been  found  that  stability  boundaries  corresponding 
to  undamped  and  slightly  damped  mathematical  structural  models 
differ  markedly.  For  example.  Loch  [l]  ,  who  considered  the  snap¬ 
ping  of  a  geometrically  perfect  sinusoidal  arch  under  a  time-step 
pressure  load,  found  a  dramatic  dynamic  influence  on  snap-through 
for  the  undamped  arch  and  certain  arch  rises;  on  the  other  hand,  his 
analysis  of  the  same  problem  assuming  a  small  amount  of  velocity- 
dependent  damping  indicates  little  or  no  difference  between  static  and 
dynamic  snap  loads  in  some  of  the  same  arch-rise  regions.  Similar 
results  were  observed  by  Huang  and  Nachbar  [2],  who  Investigated 
the  geometrically  imperfect  sinusoidal  arch  with  and  without  material 
damping  of  the  viscoelastic  type.  Further  information  is  to  be  found 
in  a  paper  by  Bieniek,  Fan,  and  Lackman  [3],  who  considered  the 
dynamic  stability  of  a  geometrically  perfect  cylindrical  shell  under  a 
spatially  uniform  lateral  pressure,  applied  as  a  step  and  ramp-step 
in  time.  There  it  was  found  that  an  undamped  analysis  predicted  a 
significant  dynamic  effect;  the  authors'  "remark",  however,  that  a 
small  amount  of  material  damping  resulted  in  the  same  stability 
boundary  for  both  the  dynamically  and  statically  applied  loads. 
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The  foregoing  investigations  indicate  the  importance  of  in¬ 
cluding  proper  damping  mechanisms  in  a  stability  analysis  involving 
transiently  applied  loads.  Further,  they  raise  some  serious  doubts 
as  to  the  validity  (without  proper  interpretation)  of  a  non-dissipative 
analysis  under  such  loading  conditions.  It  appears,  in  fact,  that  for 
a  certain  class  of  stability  problems  a  jump  condition  exists  in  the 
stability  boundary  with  respect  to  the  damping  coefficient,  say  y,  at 

y  =  0. 

In  this  paper  we  discuss  the  role  of  damping  in  stability  problems 
of  the  type  [1  -  3].  The  sinusoidal  arch  problem  of  [l]  is  employed  as 
the  vehicle  for  discussion.  This,  however,  is  a  matter  of  convenience. 
Similar  analyses  and  conclusions,  as  will  be  made  here,  can  be  made 
for  other  geometries  and  damping  mechanisms  (e.  g. ,  cylinders,  spher¬ 
ical  caps;  viscoelasticity,  etc.).  It  was  thought,  however,  that  a 
specific  example  would  better  serve  the  goal  of  this  article  than  would 
a  general  and  perhaps  undiscernible  analysis.  Thu  extension  of  the 
concepts  put  forth  here  to  other  systems  is  considered  to  be  more  or 
less  obvious. 


4 


2 


2.  Formulation 


We  consider  a  shallow,  simply  supported  sinusoidal  arch  (Fig.  1), 
which  is  governed  by  the  following  nonlinear  differential  equation  of 
motion: 

d^w  dw  a8w 

El^(w-wo)-  P  5JT+  +  P.h-jp-  *  -P(*.  t)  (2.1) 

where 

-fUL  [(&)'-(£)>  • 

Here  wo(x)  represents  the  initial  unstressed  position  of  the  arch  mid- 
surface  and  w(x)  is  the  displacement  of  said  midsurface  due  to  the 
load  p  (x,  t);  both  wq(x)  and  w(x)  are  measured  from  the  line  x  =  0 
(Fig.  1).  The  quantities  P<p  h,  E,  I,  A,  L  denote  arch  density,  thick¬ 
ness,  Young's  modulvs,  second  moment  of  area  of  the  arch  cross 
section,  cross-sectional  area,  and  length  respectively.  The  term 
multiplying  the  constant  p  is  representative  of  velocity-dependent 
damping.  The  coordinate  x  is  illustrated  in  Fig.  1 ;  t  represents 
time.  The  notation  is  that  of  [l]. 

In  addition  to  (2.  1)  we  have,  for  the  simply  supported  sinusoidal 
arch,  the  following  initial  shape  and  boundary  conditions,  respectively: 

w  =  w  sin  “■  ( w  =  const.  )  ,  (2.2) 

O  O  Li  O 

w(0,  t)  =  (0,  t)  =  w  (L,  t)  =  -gjs  (L,  t)  =  0  .  (2.3) 
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In  the  sequel  the  load,  p(x.  t),  is  assumed  to  be  of  the  form 

p(x,  t)  =  px  ^  sin  ~  ^  H  (t)  (2.4) 

where  H(t)  denotes  the  Heaviside  step  function. 

The  following  initial  conditions  are  specified  along  with 
equations  (2.  1)  to  (2.4) : 

wfx.  t^)  =  Gi(x)  ,  (x,  tQ)  =  G#(x)  .  (2.5) 

Equations  (2.  5)  can  be  physically  envisioned  as  the  result  of  a  dis¬ 
turbance  which  terminates  at  th4  time  t  =  t  . 

o 

I 

Consider  now  Figs.  2,  3.  In  [l]  Lock,  utilising  the  foregoing 
equations  of  motion,  boundary  conditions,  and  loading  function,  and 
by  virtue  of  a  numerical  integration  of  the  governing  equations  (follow¬ 
ing  their  reduction  to  two  ordinary  differential  equations  by  an  approx¬ 
imate  two-term  modal  analysis),  found  the  locus  of  the  encircled 
points  of  Fig.  2  to  be  the  snap-through  boundary1  with  a  small  value 
of  the  damping  coefficient,  y,  and  a  small  disturbance  (initial  con¬ 
ditions,  (2.5))  of  a  certain  type.  (In  Fig.  2,  e  denotes  a  nondimen- 
sional  arch  rise;  q^  is  the  nondimensional  load  amplitude  at  which 
snap-through  occurs  if  p^  is  applied  statically;  q^*  is  the  nondimen¬ 
sional  load  amplitude  leading  to  snap-through  when  p^  is  applied  as 

1  The  concept  of  dynamic  snap-through  is  defined  mathematically  in  §4. 
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*  step  in  time:  thus  q^/q  is  the  ratio  of  dynamic  to  static  snap 

( 

io&dflj.  in  particular,  note  that  little  or  no  ditterence  between  stat- 

■& 

lcally  and  dynamically  applied  loads  is  apparent  for  sufficiently  large  i 

.?X 

arch  rises.  t 

In  this  paper  we  show  that,  subject  to  sufficiently  small  dis¬ 
turbances  of  a  reasonably  general  class,  the  snap-through  boundary  for 
y  >  0  is  a  curve  in  the  q^/q^  vs.  e-plane  (Fig.  2)  which  is  defined 
by  symmetric  snap-through  up  to  the  intersection  with  the  line 

~  1*  and  q^/^  =  1  thereafter.  This  curve  is  illustrated  in 
Fig.  2  for  y  -*  0+,  and  also  for  y  corresponding  to  Lock's  1  percent 
critical  damping.  Upon  comparing  the  y  =  0+  case  with  the  y  =  0 
example  of  Lock  (Fig.  3),  one  observes  a  jump  in  the  snap-through 
boundary  at  y  =  0. 

In  contrast  to  Cl],  the  arttlysis  of  this  paper  does  not  rely  upon 
numerical  methods,  but  upon  a  rigorous  stability  analysis  incorporating 
the  full  nonlinear  equations  of  motion.  In  this  respect  the  work  contained 
herein  should  complement  [l]. 
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3-  Pr*1iminitriAi 

In  this  section  we  collect  a  few  elementary  mathematical  tools 
which  will  be  of  use  in  the  ensuing  analysis. 

Norms 

Let  the  norm,  llzll,  of  an  n- dimensional  vector  s  be  defined  as 

li*  It  =  E  !  z.|  (3.D 

1*1 

where  denote  the  elements  of  a .  Further,  let  the  norm,  ||b||  ,  of 
any  square  (n  x  n)  matrix  B  be  defined  as 

II® II  -  E  |b  |  (3.2) 

i=l 

j=l 

where  b..  are  the  elements  of  B. 
ij 

Local  Stability 

Let  t  »  0  be  a  solution  of  the  following  system  of  ordinary 
differential  equations: 

d •  (3*3) 

Here  z  is  an  n-dimensional  vector  and  g  is  an  n-dimensional  vector- 

valued  function.  We  shall  say  that  s  ■  0  is  (locally)  asymptotically 

stable  (in  the  sense  of  Liapunov)  if  1)  for  each  c  >  0  there  exists  a 

6,  depending  only  on  (  and  a  constant  r  ,  such  that  ||s(r)||  <  f  for 

o  *** 

all  r  a  r  provided  ||a(T  )||  <  fl,  and  2)  ||»(r)||-*  0  as  T  -  • . 
o  o  ~ 
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Comparison  Theorem 


The  following  comparison  theorem8  will  aid  our  discussion: 
Theorem:  Let  z  a  0  be  a  solution  of 
d  z 

—  =  A  z  +  B  (r)  z  +  f  (z ,  r )  (3.4) 

d  ^  "V  ^  IV  V  V  <v 

where  z  is  an  n-dimensional  vector,  A  is  a  constant  n  x  n  matrix 
whose  elements  depend  continuously  on  T,  and  £  denotes  a  contin¬ 
uous  n-dimensional  vector  valued  function  of  s  and  r.  (It  will  be 
assumed  that  £  is  a  nonlinear  function  of  z ,  any  linear  terms  being 
included  in  A  or  B.)  Now,  if  1)  the  trivial  solution  of 

ds 


dr 


A  7. 


(3.5) 


k 

0 

“C 


is  asymptotically  stable,  2)  B  (r)  is  impulsively  small,  i.  e. , 


f  ||B  (T)l|dT  <  M 


(3.6) 


where  M  is  a  positive  constant,  and  3)  the  function  £(js,  T)  satisfies 
the  nonlinearity  condition 


lim 

llaiho 


ll£(«.  r ) || 

W 


0 


(3.7) 


uniformly  for  r  a  0,  then  the  trivial  solution  of  (3.4)  is  asymptot¬ 
ically  stable. 


8  Theorems  of  this  type,  of  which  there  exist  many,  are  frequently  re¬ 
ferred  to  as  Poincar<- Liapunov  stability  theorems.  The  reader  is 
referred  to  Struble  [4]  and  Bellman  [5]  for  discussions. 
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Proof.  Any  solution  of  (3.4)  satisfies  the  following  integral 
equation 

*<T)  =  Y (r)  c  +  fT  Y(y  -  s)  {B(s)  s(s)  +  £(*(•).  »)}ds  (3.8) 

0 

where  Y  (r)  denotes  the  principal  matrix  solution  of  (3.  5)  and  jc  ia 
an  initial  vector.  The  hypothesis  1)  implies  that 

l|Y(r)H  *  c,  e'OT  (3.9) 

where  c^  and  a  are  positive  constants.  From  (3.  9)  and  (3.  8)  we 
obtain 

llj»(T)ll  *  cx  ||c|le"ar  +  cxJ  (  |lB(s)H  ||*(s)H  +  IU  (*(•)>  *)ll  )  <*s 

or  equivalently, 

l|s(r)|i  e“T  *  lljcU  +  ci  J  Qb(s)H  |i*(s)li  +  II £(s (s),  s)  j|  j  e®  *  ds. 

(3.10) 

Now,  the  nonlinearity  condition  (3.  7)  implies  that  there  exists  a 
positive  number  8  such  that 

li£ (a,  T)ii  *  ~  lt*ii 

l 

for  |i*||  <  fi  and  for  all  t  a  0.  Suppose  ||jcji  <  6/cx.  Then,  in  view 
of  the  continuity  of  f  ,  ||s(s)||  <  6  for  some  interval  0  s  s  *  T,  and 


from  (3.  10)  we  find 


/Till  £  r>  Ilf.  il  +  f»  !* 


on  this  interval.  Application 


f  tin  tm\\\4.  iJL.1  II  «  lm\  ll  »“*  Am 

L  .  2ClJ"~ .  - 

of  Gronwall’e  lemma  [4]  to  (3. 


/  v  in 


11)  yields 


||s(r)||  eaT  *  cjicil  «ep([  J  i|B  (s)  ||  +  2?--  ]  dr  . 

Utilising  hypothesis  2)  and  multiplying  both  sides  by  eUT  we  °k*a*n 

||.(r)ll  *  c1||cli.'‘M.-ta/2>T  .  (3.12) 

Thai,  If  ||o||  farther  eetiefiee  ||jc  ||  <  tt/c  e  1  ,  then  the  equality 

(3.  12)  will  hold  for  all  t  4  0  and  the  theorem  is  proved. 


4.  Analysis 


Basic  Equations. 

It  will  be  convenient  ,  as  in  [l  3,  to  nondimens ionalize  equations 
(2.  1)  to  (2.  5)  by  introducing  the  new  variables 


(4.  1) 


where 


.  nx 

4  =  L  ' 

r  = 

“o* 

:3 

n 

XI* 

w 

•a 

o 

II 

w 
_ 0 

k 

r8 

/  El  1 

Nl/# 

and  k  is  the  cross-sectional  radius  of  gyration.  With  use  of  (4.  1), 
equations  (2.  1)  to  (2.  5)  reduce  to 

5?- -  v  ■  &  $ !”[(!? )a-  (^)”]d< + yl? +  & 1  -,(t' T) 


where 


q(C,  T)  =  ^  sin  i  H(r) 


JL(k 

%.  ~  Elk  \  ir 


p  h  W 


(4.2) 


(*)  • 


(4.3) 


In  addition,  we  have 


t?q  =  e  sin  $ 


e  =  w/k  , 
o 


-3  - 

n<0.  t)  =  §”?<o,t)  =  n<*.  T)  =  ^r(^r)  =  0  , 

T ?<{.  0)  -  gjU)  ,  0)  =  g8(0  , 
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(4.4) 

(4.5) 

(4.6) 


where  =  G^/k,  i  =  1,2. 

An  exact  solution  of  equation  (4.  2),  satisfying  the  boundary 
conditions  (4.5)  termwise,  is 


T?(£.  r)  =  e  sin  £  +  £  a  (r)  sin  n$ 


(4.7) 


where  the  coefficients  aR  (r)  are  solutions  of  the  following  system 
of  ordinary  differential  equations: 


»’•»  ) <*»  +  •>  ' 

n=l 


(4.8) 


+  y  —t~  +  m4a  +  j  (  2ea  +  £  naa  8  )  ~  0  • 

dr  m  4  \  i  ,  n/  m 


m  —  2, 3, . . . ,  N 


To  avoid  the  question  of  convergence,  we  shall  assume  that 
the  number  N  in  (4.  7)  and  (4.  8)  is  finite.  This,  of  course,  restricts 
the  class  of  disturbances  that  the  arch  may  experience.  In  particular, 
(4.  7)  constitutes  an  exact  solution  only  if  g^  and  ga  can  be  expressed 


in  the  form 


N  da 


*i  ({)  =  E  an(rQ)  sin  n$  ,  g^)  =  £  (To)  gin  ^  .  (4.  9) 


The  initial  conditions  to  (4.  8)  are  the  Fourier  coefficients 


W  =  v  j0  *.<«> 


"*»<%>  1,' 

— dr - ”  «, («)•!■>  (4.10) 


.**!■*■« 


From  a  physical  point  of  view,  the  above  finite  N  restriction  is  not 

considered  to  be  overly  important. 

Symmetric  Motion  and  Symmetric  Snap-through. 

Under  the  assumption  that  g^(£)  5  0,  i  =  1,  2,  the  response  of 

the  arch  is  purely  symmetric,  i.  e. ,  a  *  0,  n  =  2,  3,  .  .  . ,  N  and  a 

n  l 

is  governed  by 


d*ax  dax 

+  yd7  +  ax  +  4  (2eV  ata>  <V  e)  s  . 


*x(0) 


dal(0) 

~dT 


0 


(4. 11) 


Now,  let  a^  denote  the  equilibrium  position(s)  the  arch 
assumes  if  is  applied  statically.  Then,  a^  satisfies 

a  +  y  (2ea  +  a  a )  (a  +  e)  +  a  -  0  .  (4.  12) 

e  4  e  i  e  i 

For  <  (51)g»  and  for  e  >  2,  there  exist  three  real  roots  of  the 
cubic  (4.  12)  which,  in  turn,  define  three  equilibrium  positions  of  the 
arch.  For  qi  >  (IJ^  )g,  there  exists  only  one  real  root  and  hence  only 
one  equilibrium  position  (the  inside-out  position  of  the  arch).  The 
number  (q^  )g  is  the  classical  static,  symmetric,  snap-through-load. 

Let  us  assume  that  ql  <  and  let  us  introduce  the  trans¬ 

formation 

ai  (r)  =  ai  "  V  <T>  (4.  13) 

where  a  =  max  a  ,  i.  e.  ,  the  unsnapped,  symmetric  equilibrium 
•  © 


position. 
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By  virtue  of  (4.  11),  a^'  '  (r)  satisfies 


dV(0) 


+  3(a(0);a)  =  0 


(4.  14a) 


a^0,(0)  =  ax  , 


(0)  =  0 


(4.  14b) 


where 


5  *  •1(0,+  4  +  3V  +  ^ J  -  *i°)a[3<I,+  «'J  +  *,(0>a}  ' 

(4.  14c) 

In  terms  of  the  new  variable  a^°\  the  equilibrium  positions  a# 
are  the  critical  points  of  (4. 14),  i.  e. ,  da^^/dr  =  0,  *  (a  ^)j, 

I*  1,  2,3,  where  the  (aj0^  are  defined  by  3  <ai<0);  I)  =  0  and 
ordered  such  that  (a(0)),  <  (a(0))a  <  (a(0))  .  The  point  (a,(0))  is 
the  origin  a^®^  *  0  and  corresponds  to  the  unsnapped  configuration 
of  the  arch;  the  point  (a  ^)  is  the  inverted  arch  position. 

Consider  next  the  phase  plane;  da^/dr  vs.  a  In  the 
phase  plane  it  can  be  shown  that  (a^*^  is  a  center  point,  (a^^  is 
a  saddle  point,  and  (a^)  is  another  center  point.  If  q^(i.e.  ,  the  initial 
data  (4.  14b))  is  such  that  the  trajectories  in  the  phase  piano  encompass 
the  saddle  point  (»1<°))a ,  we  shall  say  (.ee  [8])  that  dynamic,  sym¬ 
metric,  snap-through  has  occurred.  On  the  other  hand,  if  all 
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!•  -  -  -  -  .  -1  .  .1  _1.  .  J  ll  ■  •  1  -1  -  ..  .i  ......  .1  ^  .  -*■  1*  -  -  -■*  —  ' 

biujvvvv  *  it  a  ai  t  tiuotu  duuut.  u  i  it;  i/xigui  aau  uu  uwt  caiciiu  *\j  v*  wv  y 

the  saddle  point,  we  shall  say  that  dynamic,  symmetric,  snap-through 
i:>i3  not  occurred.  The  value  of  q  leading  to  a  trajectory  which  inter¬ 
sects  <~he  saddle  point  represents  the  transition  between  snapped  and 
unsnapped  ar.-hes;  we  shall  denote  this  value  as  (c^*)  . 

In  addition  to  the  assumption  that  <  (q^)^,  let  us  ai,ume  *1^ 
is  leas  than  that  required  to  cause  dynamic,  symmetric,  snap-through,  i.  «. , 
qi<  (q*)  .  Then,  for  y>  0,  the  trajectories  in  the  phase  plane  spiral  toward  the 
origin.  Thus,  for  a  given  qx  and  y,  and  for  any  c  >  0,  there  exists  a  num¬ 
ber  T  (f  ;  q)  such  that  |a^*(T  )|  +  jda^/df(T  )|  <  C  for  all 
T  s  T  .  Further,  since  the  nonlinear  portion,  n3,  of  the  function  3 
satisfies  the  nonlinearity  condition 

lim  In  3  | 

|.'°'l-0  I.*01!  =  ’ 

a  well-known  comparison  theorem  [4]  indicates  that 


l*1(0)(T)|  *  Pe 


-(y/2  )T 


,  p  =  const. 


(4.15) 


for  sufficiently  small  c,  or  for  sufficiently  large  T  .  From  (4.  15), 
we  therefore  conclude  that  a^(r)  is  impulsively  small  as  T  •  , 


03 

[  |a  ^  (r ) I  dr  <  M  =  const. 

V1 


(4.  16) 


We  shall  make  use  of  this  fact  later. 


Definition  of  Dynamic  Snap-through 

The  concept  of  dynamic  snap-through  has  been  mentioned  in 
connection  with  the  purely  symmetric  motion.  In  order  to  avoid  a 
vague  discussion,  we  must  attempt  to  pin  down  this  illusive  butterfly 
for  a  .symmetric  as  well  as  symmetric  motions.  The  following  defini¬ 
tion  will  suffice  for  purposes  of  the  present  analysis: 

Definition.3  Let  us  perturb  the  symmetric  motion  according  to 


\(T)  =  \  (T)  +  a^*  (t ) 


»m(r)  =  am*(T)  ,  m=2,  i . N. 


(4.17) 


Next,  let  us  assume  that  qi  is  given.  We  shall  say  that  snap-through, 

under  load  q^,  has  not  occurred  if  there  exists  a  6  >  0  such  that 
N 

E  )|  +  Ida  */dr(T  )j  }  <  8  implies  max  [s  _a  (r)]<  (a*0)) 

n=l  (Xt<®  11  1  J 

for  any  constant  Tq  (recall  that  (a  ^)^  is  the  inside-out  equilibrium 

position).  Conversely,  we  shall  say  that  snap-through  has  occurred  if 

there  exists  at  least  one  set  of  initial  data  such  that 

max  [a  -  a  (T)  3  a  (a  ^)  for  arbitrary  small  6. 

0<T<®  11  13 

Local  Stability  of  Symmetric  Motion. 

Let  us  now  investigate  the  local  stability  of  the  symmetric  motion 


3  This  definition  is  compatible  with  the  snap-through  concepts  employed 
in  both  [  1  ]  and  [  2] . 


\lr 


IU  J,  \ 

aj  w  )  with  respect  to  all  possible  symmetric  and  asymmetric  motions 
resulting  from  the  disturbance  (4.  10).  To  accomplish  this  we  perturb 
a  ^(T)  according  to  (4.  17).  Upon  substitution  of  (4.  17)  into  (4.  8)  one 
finds  that  the  perturbations  a^  satisfy  the  following  system  of 
equations: 


i  da 

s'  *  ’"5T  +  \*  [‘  +  4  (6“,  +  3V  *  )j 


3a*a(0)(T)  3a  ** 

+  — - [»l(0)(T)  -  2(1  +  e)j  +  -j—  (I  +  e) 


‘l*  +  *>  N 

+  -j-  +  - j -  E  n8a  **  - 


*1 0)  (T)  N  ,  a 
-  £  n8  a  ** 

4  n=2 


3a* 

-T-  •> 


(4.  18a) 


*i*  N 

+  -T-  £  naa  *a  .  0  , 

*  *  n 


d8a  * 
m 


da  *  r  a,m8  -I  *  a  0  (T)m8  . 

v]  *  —  [s,0,«  - 


-  2(1 


:1+  *>] 


(4.  18b) 


2  m  l 


4  m 


+  ^  a  *  g  n8  a  * 
4m.  m 
n=2 


0  ,  m  :  2 .  N, 
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With  the  substitutions 


(It, 3  Hi 

dT 


'i  “  (i/2) 


i  =  1,3,5,...,  2N- 1 


,  i  =  2,  4,  6,  ....  2N 


(4.  19) 


equations  (4.  18)  can  be  reduced  to  the  following  2N  system  of  first 


order  nonlinear  differential  equation*: 


+  5<r>£,  +  i(*>  r> 

dT 


(4.  20) 


Here  «  denotes  a  2N- dimensional  vector  whose  elements,  *.,  are 
^  1 

defined  by  (4.  1 9);  A  is  a  2N  X  2N  constant  matrix,  the  form  of 
which  is  not  pertinent  to  our  discussion;  B  (T)  is  a  2N  x  2N  matrix 


having  the  form 


B  (T)  *  *}0)  (T)  [V0)  (r)  -  2(ax  -  e)J  C  (4.  21) 


where  C  is  a  2N  x  2N  constant  matrix.  The  2N-dimen*ional 
vector  f  in  (4.  20)  consists  of  autonomous  (say  f  ^)  and  nonauton- 
omous  (say  )  parts.  Writing  J  =  f^  +  we  have: 


3z  J  z‘3  ( a  +  e  +  z  )  m 

<i  ♦  .)  ♦  -f  «  — - — r — -  s 

l  4  i  4  4 

n=2 

-  f.(l5  =  0  ,  i  =  2,  4,  6,  .  .  .  ,  2N  , 


_a  _  3 
L  n  z 

an 


(4.  22a) 


a  a  a  M 

,h\  m  .  —  .  .  m  a  ,  rn  _  _a  _ 

f.  ’  =  ~T”  a  a  ( a  +  e)  +  - t~  a  *  +  “T"  a  L  n  a 

i  2  a  am  '  i  4  a  am  4  am  r_2 


.a  _  a 
am 


i  =  3,  5,  7.  .  .  .  ,  2N-1;  m  =  <i  +  l)/2  . 


fj3) 


*i°\t) 


K 


N 

I 

n=2 


na  a 


al 

an  J 


f(a) 

i 


=  0  ,  i  =  even 


(4.  22b) 


f(a)  a  IE.  z  B  a(0,(T)  ,  i  =  3,5,7,...,  2N- 1  ; 

i  2  am  a  l 

m  =  (i  +  1 ) /2  . 


Let  us  now  assume  that  ^  <  (q^,  ^  <  t^*)  .  Then,  by 
virtue  of  (4.  1 6),  we  conclude  that 

J"  ||b(t)||  dr  <  M  =  const.  (4.23) 

0 

Further,  noting  that  a  ^(t)  is  a  bounded  function  of  r,  and  (from 
(4.  22))  that  the  elements  of  £  are  polynomials  in  e^  with  no  linear 
terms,  it  is  evident  that 


Urn  ||f  ||  /  ||a  ||  =  0 

llsil-0  ~  ~ 


(4.  24) 


Applying  now  the  comparison  theorem  of  Section  3,  we  conclude  that 
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z  e  o  (the  symmetric  motion)  is  asymptotically  stable  provided  the 
origin  of 

d  *, 

— —  =  A  z  (4.  25a) 

dr  ~ 

is  asymptotically  stable. 

In  terms  of  the  original  system  of  N  second-order  equations, 

(4.  25a)  can  be  written 


d3a  * 

,  # 
da 

l 

_______  i 

*r,  1 

ai  L1  +  4 

/  L  „  T-  i 

<5  ‘T"  3  i  d  _« 

y 

dr 

I  o  e  a  t 

V  l 

j  a  t  c  c 
l 

d8a  * 
m 

da  * 

_ m  + 

a  *  [“  a 
m  a  I  m 

a 

+  jl  (2e 

+  i)l  = 

dr8  + 

y 

dr 

m  L 

1  j 

(4.  25b) 


m  =  2,  3,  . . . ,  N. 


One  recognizes  (4.  25b)  as  the  classical  perturbation  equations  for  the 
statically  loaded  arch  [l,  6,  7].  Traditionally,  the  critical  static  snap- 
through  load,  ,  is  based  on  such  equations.  In  particular,  if  y  *  0, 
corresponds  to  the  transition  between  bounded  and  unbounded 
solutions  of  (4.25).  If  y  >  0,  and  <  q^ ,  all  solutions  of  (4.25) 
decay  exponentially  to  zero.  Thus,  by  definition  of  ,  if  q^  <  q^ 
and  y  >  0,  the  origin  of  (4.  25)  is  asymptotically  stable.  We  thus 
conclude  that  the  origin  of  (4.  20)  is  also  asymptotically  stable  if 

%.  <  V  qi<  (v.’  qi<  {t*iV 
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Local  Stability  and  Snap- through. 

Consider  now  the  foregoing  result  in  terms  of  snap-through. 
Local  asymptotic  stability  implies  that,  under  sufficiently  small  distur¬ 
bances,  all  motions  of  the  arch  can  be  maintained  within  an  arbitrarily 
small  neighborhood  of  the  symmetric  motion.  To  be  more  precise,  we 
have  the  following:  if  y  >  0  and  qj  satisfies  i)  q^  <  (q^,  ii) 

\  <  (qi*)8.  iii)  ,  then  for  any  €  >  0,  there  exists  a  6  >  0 

(depending  in  general  on  Tq,  y,  q^  and  the  arch  geometry)  such  that 

||s(t)H  <  c  for  all  r  2  r  provided  ||*(r  )i|  <  5.  Thus,  under  such 

o  o 

conditions  asymmetric  snap-through,  as  defined  previously,  is  not 
possible.  One  therefore  concludes  that,  if  y  >  0  and  q^  lies  in  the 
region  defined  by  i) ,  ii),  iii)  above,  then  for  sufficiently  small  distur¬ 
bances  snap-through  will  not  occur.  On  the  other  hand,  if  >  (qi*)|i* 
it  is  evident  that  symmetric  snap-through  will  take  place.  Further,  the 
analysis  of  Huang  and  Nachbar  [2]  implies  there  exists  a  set  of  initial 
conditions  (resulting  in  an  excitation  of  only  two  modes,  i.  e.  ,  for  T  =  0, 


(0) 


(0) 


am  =  dam  /dT  =  0  for  m  >  2)  such  that  snap-through  takes  place 
for  an  arbitrarily  small  initial  disturbance  if  qi  >  q^  and  y  >  0. 
Therefore  the  stability  boundary  for  infinitesimal  disturbances  must  be 
a  curve  consisting  of  dynamic,  symmetric  snap-through  up  to  its  inter¬ 
section  with  <^*7^  -  1.0  and  qi  =  1.0  beyond.  Finally,  it  is  evident 
that  this  curve  also  constitutes  an  upper  bound  on  stability  for  distur¬ 
bances  of  finite  amplitude. 
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The  foregoing  curve  is  illustrated  in  Fig.  2  as  y  —  0+,  and  for 


y  corresponding  to  \°!z  critical  damping  in  the  symmetric  mod* 


(Lock's  case).  The  slight  disagreement  with  Lock's  numerical  results  is 


attributed  to  the  siae  of  the  initial  conditions  employed  in  that  analysis, 


i.  e. ,  presumably  the  smaller  the  initial  data,  the  better  the  agreement 


would  be.  The  symmetric  snap- through  portion  of  the  curves  exhibited 


in  Fig.  2  was  obtained  from  data  supplied  by  Lock  for  y  f  0  (numerical 


integration  of  the  differential  equation  (4.  11)  is  necessary  if  y  f  0. ) 


An  exact  relation  for  (q^*)^  vs.  e  can  be  obtained  for  y  =  0+,  e.  g. 


from  [2].  In  the  latter  case  a  combination  of  Eqs.  (32),  (35)  and  (48) 


of  [2]  yield  the  following  relations  for  the  snap- through  boundary 


(dotted  lines  of  Fig.  2)  when  y  =  0+: 


(i)  e  <  /IT  =  4.68 


(1/27)  [e<9  +  y  )  *  (£  -  3)  (•’  -  ] 


2  +  f  /T  <4  ‘  »V> 


,  (4. 26a) 


(ii)  4.68  s  e  *  6.38 


(1/27)  [e(9  +  y)  +  ( y  -  3)  (ea 


-  6)V8  ] 


f  +  S(£  -  4)V» 


,  (4. 26b) 


(ill)  e  a  6.38 


%*'S  - 


(4.  26c) 


5.  Concluding  Femarks 


As  an  example  ol  a  structure  exhibiting  a  particular  sensitivity 
to  damping,  the  snapping  of  a  simp’y -supported,  geometrically  perfect 
sinusoidal  arch,  subjected  to  a  sinusoidal  spatial  pressure  distribution 
applied  as  a  step  in  time,  was  considered.  It  was  found  that,  in  the 
presence  of  velocity-dependent  damping  of  any  magnitude  (>0),  and  under 
sufficiently  small  disturbances  (taking  the  form  of  initial  arch  displace¬ 
ments  and  velocities),  there  exists  no  difference  between  dynamic  and 
static  snap-through  pressures  for  sufficiently  large  arch  rises.  For  in¬ 
finitesimal  disturbances  the  snap-through  boundary  is  a  curve  in  a  load 
vs.  arch  rise  plane  defined  by  symmetric  snap  up  to  its  intersection  with 
the  static  snap-load,  and  the  static  value  thereafter.  For  finite  distur¬ 
bances  this  represents  an  upper  bound  on  this  load.  Said  curve  is  illus¬ 
trated  in  Fig.  2  for  vanishingly  small  damping.  Upon  comparison  with 
the  purely  elastic  case  of  Lock  [l],(Fig.  3),  one  observes  a  jump  in  the 
snap-boundary  with  respect  to  damping. 

A  similar  analysis  as  has  been  conducted  here  can  be  applied  to 
other  structures.  In  particular,  an  analysis  of  geometrically  perfect 
rings  and  cylinders  under  spatially  uniform  loading  applied  as  a  step  in 
time  will  yield  the  result  that,  under  sufficiently  small  disturbances, 
there  exists  no  difference  between  statically  and  dynamically  applied 
buckling  loads  if  damping  is  present. 

Finally,  we  close  on  an  interesting  note.  Suppose  the  foregoing 
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load  were  applied  as  an  impulse  rather  than  a  step;  then  a  *  0.  In 
this  case  (4.  20)  is  asymptotically  stable  for  all  .  Therefore,  one 
excludes  that,  in  the  presence  of  damping,  and  under  sufficiently  small 
disturbances,  only  symmetric  snapping  is  possible.  In  the  case  of  the 
ring  and  cylinder  mentioned  previously,  no  dynamic  buckling  whatsoever 
would  result.  The  physical  explanation  of  this  is  as  follows;  in  the 
presence  of  damping  a  sufficiently  small  disturbance  is  damped  out 
before  it  can  initiate  a  "mode  conversion"  (by  virtue  of  parametric 
resonance)  which  results  in  a  buckling  phenomena.  The  foregoing 
examples  indicate  that  the  magnitude  of  both  disturbances  and  damping 
are  important  for  near  perfect  structures  under  certain  transient  loadings. 
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